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The form factors of the weak Bg transitions to the ground state and orbitally 

excited strange mesons are calculated in the framework of the QCD-motivated rela- 

fr^ \ tivistic quark model based on the quasipotential approach. These form factors are 

^ expressed through the overlap integrals of meson wave functions found in their mass 

^vg ' spectrum evaluations. The momentum dependence of the form factors is determined 

^ ■ in the whole accessible kinematical range without any additional assumptions and 

extrapolations. Relativistic effects, including the wave function transformation from 

the rest to a moving reference frame as well as contributions of the intermediate 

negative-energy states, are consistently taken into account. The calculated form 

factors are used for the evaluation of the charmless semileptonic decay rates and 

' (— i '' two-body nonleptonic Bs decays in the factorization approximation. The obtained 

Qh' results are confronted with previous predictions and available experimental data. 
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I. INTRODUCTION 
> 

in 

\f^ • In last years the study of the properties of B^ mesons and especially its weak decays 

^ . attracts attention of both theorists and experimenters. Indeed the investigation of such 

decays is important for the independent determination of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements, studying the CP-violation, testing the standard model and a "new 
cn ' physics" models. New precise data are coming from the experiments on Large Hadron 

Collider (LHC) which significantly extended the number of the observed Bg decay channels 

^^ \ In the recent paper [Sj] we considered the weak Bs decays to the ground and excited 

H ' states of charmed mesons. Such decays are the dominant decay channels of the Bg meson. 

- - - We calculated the weak decay form factors in the framework of the relativistic quark model 

with the QCD-motivated interquark potential. Good agreement of our predictions for the 
semileptonic and nonleptonic decay branching fractions with experimental data has been 
found. Here we extend our analysis to the consideration of the CKM suppressed Bg decays 
to strange mesons. First we calculate the corresponding weak decay form factors. The 
characteristic features of the heavy-to-light Bs ^ K transitions are the presence of the light 
K meson in the final state and a very broad accessible kinematical region. Therefore it is 
important to consistently take into account the relativistic effects, which give substantial 
contributions, and to reliably determine the momentum dependence of the form factors. 
It is necessary to point out that most of the available theoretical approaches permit to 
evaluate the form factors either at some fixed point (of zero or maximum recoil of a final 
light meson) or determine the momentum dependence of the form factors in a restricted 
kinematical range. As a result they require ad hoc assumptions about the q^ dependence of 
form factors or their extrapolations, which also rely on specific models. The advantage of 



TABLE I: Masses of the ground state and first orbitally excited strange mesons calculates in our 
model (in MeV). 
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our model consists in its ability to explicitly determine the corresponding decay form factors 
as the overlap integrals of meson wave functions in the whole kinematical range without any 
additional assumptions and extrapolations. 

In our calculations of the decay form factors we use the wave functions of the Bg and K 
mesons obtained previously in their mass spectra evaluations [J, l5j. In Table [I] we compare 
our predictions for the masses of the ground state and lowest orbitally excited strange mesons 
[Sj with available experimental data J2|. We find satisfactory agreement of our predictions 
with data. Note that in our mass spectrum calculations all relativistic effects, including the 
spin-dependent and spin-independent contributions to the potential, were treated nonper- 
turbatively in v^ !<? . From Table [T] we see that our model predicts that all orbitally excited 
Kj mesons have masses heavier than 1 GeV. The scalar i^o(800) (or /^) meson is predicted 



in our model to be a scalar tetraquark [6|, while some other theoretical approaches assume 
it to be the scalar quark- ant iquark (1^-Po) state. It is clear that the form factors of the weak 
Bs — !■ i^Q transition are significantly different if the ii'Q(1430) meson is the I^Pq or 2^Po 
state. Therefore the resulting rates of the semileptonic and nonleptonic B^. decays to the 
i^o(1430) strongly vary depending on its structure and quantum numbers. Thus the study 
of the weak B^ decays to the K^ meson can help to reveal the nature of light scalar mesons. 
In the following we denote the iTgll^SO) meson by K^. 

The calculated form factors are used for evaluation of the rates of the charmless semilep- 
tonic decays both to the ground state and orbitally excited strange mesons. The two-body 
tree-dominated nonleptonic B^ decays are considered within the factorization approxima- 
tion. The charmless B^ decay rates are evaluated and compared with previous calculations 
and available experimental data. 



II. RELATIVISTIC QUARK MODEL 



The employed relativistic quark model is based on the quasipotential approach in quan- 
tum chromodynamics (QCD). Hadrons are considered as the bound states of constituent 
quarks which are described by the single-time wave functions satisfying the three-dimensional 
relativistically invariant Schrodinger-like equation with the QCD-motivated interquark po- 
tential p] 
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where the relativistic reduced mass is 

^"^ ~ 4M^ ' ^^' 

M is the meson mass, mi 2 are the quark masses, and p is their relative momentum. In 
the center of mass system the relative momentum squared on mass shell 6^(M) is expressed 
through the meson and quark masses: 

_ [M^-{m^ + m2f][M^-{rm-m^f] 

The kernel of this equation is the interquark quasipotential V{p, q; M) which consists of the 
perturbative one-gluon exchange and the nonperturbative confining parts [7| 

V{p, q; M) = u,{p)u2{-p)Vip, q; M)ui{q)u2{-q), (4) 

with 

V(p, q; M) = ^asD,Ml^l2 + Ko„f (k)r^(k)r2;^(k) + KL(k), k = p - q, 

where a^ is the QCD coupling constant, Dfj^u is the gluon propagator in the Coulomb gauge, 
and 7^ and u{p) are the Dirac matrices and spinors, respectively. The Lorentz structure 
of the confining part includes the scalar and vector linearly rising interactions which in the 
nonrelativistic limit reduce to 

Konf (r) = Konf (0 + Konf (0 = Ar + B , (5) 

with 

Konf(^) = (1 - ^){Ar + B), Konf(0 = e{Ar + B), (6) 

where e is the mixing coefficient. Its value e = —1 has been obtained from the consideration 
of the heavy quark expansion for the semileptonic B ^ D decays |8| and charmonium 
radiative decays [7(]. 

The long-range vector vertex 

r,(k) = 7, + ^cr^uk^ (7) 

contains the Pauli term with anomalous chromomagnetic quark moment k. The value k, = 

— 1, fixed in our model from the analysis of the fine splitting of heavy quarkonia ^Pj- states 
7 1 and the heavy quark expansion for semileptonic decays of heavy mesons [8] and baryons 
9|, enables vanishing of the spin-dependent chromomagnetic interaction, proportional to 

(1 + k), in accord with the flux tube model. 

Other parameters of our model were determined from the previous analysis of meson 

spectroscopy [7\. The constituent quark masses are rub = 4.88 GeV, rric = 1.55 GeV, 

rUs = 0.5 GeV, m„,rf = 0.33 GeV and the parameters of the linear potential are A = 0.18 

GeV^ and B = -0.30 GeV. 

For the consideration of the meson weak decays it is necessary to calculate the matrix 

element of the weak current between meson states. In the quasipotential approach such 
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FIG. 1: Leading order vertex function r(^^(p, q) contributing to the current matrix element ([8|). 

matrix element between a Bg meson with mass Mb^ and momentum pbs and a final K 
meson with mass Mk and momentum px is given by [lO| 



where r^(p, q) is the two-particle vertex function and ^ mpm^v) are the meson (M = Bg, K) 
wave functions projected onto the positive energy states of quarks and boosted to the moving 
reference frame with momentum pj\/, and p, q are relative quark momenta. 

It is convenient to carry calculations in the Bg meson rest frame (p^^ = 0). Then the 
final meson is moving with the recoil momentum A. The wave function of the moving meson 
^ K/s. is connected with the wave function in the rest frame ^ ko = "^k by the transformation 

i3 

^K^ip) = Dl/\RYjDl/\RYj^Ko{p), (9) 

where R^ is the Wigner rotation, La is the Lorentz boost from the meson rest frame to a 
moving one and D^^'^{R) is the spin rotation matrix. 

The wave function of a final '^^^^Kj meson at rest is given by 

^Kiv) ^ ^/i;^.^,(p) = y'"-'^ ^^.+ii.,(p), (10) 

where J and M. are the total meson angular momentum and its projection, L is the orbital 
momentum, while S* = 0, 1 is the total spin. '?/'2s+iXj(p) is the radial part of the wave 
function. The spin-angular momentum part y^^^^-'^ is defined by 

0-10-2 -^ ^ 

(11) 

where {jimi, J2^T^2\J -M) are the Clebsch-Gordan coefficients, Y[^ are the spherical har- 
monics, and x{(^) (where a = ±1/2) are the spin wave functions. 

The explicit expression for the vertex function r^(p, q) can be found in Ref. |3|]. It con- 
tains contributions both from the leading order spectator diagram (see Fig. [1]) and from 
subleading order diagrams accounting for the contributions of the negative-energy interme- 
diate states. The leading order contribution 

rW(p, q) = u,{p,h,{l - 7')Mb(g,)(27r)3,5(p, - q,) (12) 

contains the 6 function which allows us to take one of the integrals in the matrix element ([8]). 
Calculation of the subleading order contribution is more complicated due to the dependence 



on the relative momentum in the energies of the initial heavy and final light quarks. For 
the energy of the heavy quark we use heavy quark expansion. For the light quark such 
expansion is not applicable. However, the final light K meson has a large (compared to its 
mass) recoil momentum (lAmaxI = {^b ~ ^k)/{'^-^Bs) ~ 2.6 GeV) almost in the whole 
kinematical range except the small region near q^ = q^^^^ (|A| = 0). This also means 
that the recoil momentum of the final meson is large with respect to the mean relative 
quark momentum |p| in the meson (~ 0.5 GeV). Thus one can neglect |p| compared to 
|A| in the light quark energies eq{p + A) = 



m„ 



IP 



replacing it by eq(A) 



ml + A^ in expressions for the subleading contribution. Such replacement removes the 
relative momentum dependence in the energies of quarks and thus permits to perform one 
of the integrations in the subleading contribution using the quasipotential equation. Since 
the subleading contributions are suppressed the uncertainty introduced by such procedure is 
small. As the result the weak decay matrix element is expressed through the usual overlap 
integral of initial and final meson wave functions and its momentum dependence can be 
determined in the whole accessible kinematical range without additional assumptions. 



III. FORM FACTORS OF THE WEAK TRANSITIONS OF B^ TO K MESONS 

The matrix elements of the vector, axial vector and tensor weak currents between Bs and 
K^*^ meson states are parametrized by the following set of form factors 
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q = pb, —Pk(*) 5 ^^(^ Mb k(*) are the masses of the B meson and final K^*^ meson, respectively; 
while e^ is the polarization vector of the final vector K* meson. 



TABLE II: Calculated form factors of weak Bg — )• K^*' transitions. Form factors f+{q^), /t(9^) 
y(g2), ^o(g'), ri((?2) are fitted by Eq. ([lOD, and form factors Uq^), Ai{q^), ^2(9'), ^2(^2), T^{q^) 
are fitted by Eq. ^. 
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At the maximum recoil point (g^ = 0) these form factors satisfy the following conditions: 

/+(0) = /o(0), 
, , , Mb + Mk* ... Mb - Mk* . . . 

Ti(0) = T2(0). 

Comparing these decompositions with the results of the calculations of the weak current 
matrix element in our model, as described in the previous section, we determine the form 
factors in the whole accessible kinematical range through the overlap integrals of the meson 
wave functions. The explicit expressions are given in Refs. [ll|, |l2|. For the numerical 



evaluations of the corresponding overlap integrals we use the quasipotential wave functions 
of Bg and i^*^*-* mesons obtained in their mass spectra calculations j^ |5|. 

We find that the weak Bg — )■ K^*^ transition form factors can be approximated with good 



accuracy by the following expressions [l3|, |l4 1 
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where M = Mb* for the form factors /+(g^), /r(g^), 1^(g^),Ti(g2) and M = Mb for the 
form factor Aolq"^). The obtained values F{0) and (Ti^2 are given in Table HTl The difference 
between fitted and calculated form factors is less than 1%. We can roughly estimate the total 
uncertainty of the form factors within our model to be less than 5%. It mainly originates 
from the subleading contributions to the decay matrix elements in the region of small recoils. 
These form factors are plotted in Figs. [2] and [31 

In Table IIIII we compare our predictions for the form factors of weak Bs decays at max- 



imum recoil q = with results of other calculations |13l . |15|-|21||. Light-cone sum rules. 
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FIG. 3: Form factors of the weak Bg — )• K* transition. 



including one-loop radiative corrections to twist-2 and twist-3 contributions, and leading 
order twist-4 corrections are used in Ref. jl5|. Perturbative QCD approach is applied in 
Refs. 16|, llTi ■ Calculations based on quark model and relativistic dispersion approach are 
given in Ref. 13|, while consideration in Ref. |l8| is performed in the light cone quark model 
utilizing the soft collinear effective theory. The authors of Re f. |l9l| employ light-cone sum 
rules in the framework of heavy quark effective theory. In Ref. [20| the weak transition form 
factors are evaluated in the six-quark effective Hamiltonian approach. Perturbative QCD 
factorization approach with the inclusion of the next-to-leading-order corrections is used 
in Ref. 2l|. We find a reasonable agreement between the values of the weak Bs — ^ K^*'^ 
transition form factors at g^ = calculated in significantly different approaches. 



IV. FORM FACTORS OF WEAK TRANSITIONS OF Bs MESONS TO 
ORBITALLY EXCITED K^^ MESONS 

Now we apply the same approach for the calculation of the form factors of the weak Bs 
decays to orbitally excited Kj mesons. The matrix elements of the weak current J^ = 
^7^(1 — 75)6 for Bs decays to orbitally excited P-wave Kj mesons can be parametrized by 



TABLE III: Comparison of theoretical predictions for the form factors of weak Bg — >■ K^*' transitions at maximum recoil point q"^ = 0. 
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the following set of invariant form factors 

{K*,ipK,)\urb\B,ipBj) = 0, (22) 

{K*{pK,)\ur^,b\B,{psJ) = r^iq') (p^, +P^J + r_(g') {pk ' Pk) ' (23) 
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(ir*(p;,J|M7'^75fe|5s(j>Bj) = (Mb, + MK,)tAAqn^ ^ 

+ [tM{q'VB. + tA,{qVKK^^^, (27) 

where g = pb^ — Pkj, ^Kj are P-wave K meson masses, e^ and e^'^ are the polarization 
vector and tensor of the vector Ki = i^i(1270) and tensor K2 mesons, respectively. The 
matrix elements of the weak current for Bs decays to the axial vector ii'i(1400) meson are 
obtained from Eqs. (!24l) . fl25l) by the replacement of the set of form factors hi{q^) by gi{q^) 

The P-wave K meson states with J = L = 1 are the mixtures of spin-triplet (^Pi) and 
spin-singlet {}Pi) states: 

\Ki{l27Q)) = \KCPi))cosip+\KCPi))sm(p, 

|Xi(1400)) = -|ir(iPi))sin¥P+|ir(^Pi))cos(^, (28) 

where y) is a mixing angle. Such mixing occurs due to the nondiagonal spin-orbit and tensor 
terms in the spin-dependent part of the relativistic quasipotential. The masses of physical 
states are obtained by diagonalizing the mixing terms. The found value of the mixing angle, 
if = 43.8° [5], implies that physical Ki mesons are nearly equal mixtures of the spin-singlet 
K{^Pi) and spin-triplet K{^Pi) states in accord with the experimental data [2!]. 

We calculate the matrix elements of the weak current between the initial P» meson and 
final orbitally excited P-wave Kj meson using the procedure described above and compare 
the result with the invariant decomposition f l22|) -f l27|) . In this way we obtain expressions 
for the decay form factors in terms of the overlap integrals of the initial and final meson 
wave functions which are valid in the whole accessible kinematical range. Note that these 
form factors account for the relativistic transformations of the meson wave functions from 
the rest to a moving reference frame (Q and relativistic contributions from the intermediate 
negative energy states. The explicit expressions for these form factors can be obtained from 



the corresponding formulas in Appendix of Ref. [llj with obvious replacements. 

For the numerical evaluations we again use the meson wave functions obtained in their 
mass spectra calculations j^, |5[. The calculated values of Bs to orbitally excited Kj transi- 
tion form factors at maximum and zero recoil are give in Table IIVI The momentum depen- 
dence of these form factors is shown in Fig. m We can estimate the errors in the calculated 
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TABLE IV: Form factors of the weak Bs decays to the P-wave K j mesons at q^ = and 
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FIG. 4: Form factors of the Bs decays to the P-wave K K mesons. 
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form factors to be less than 10%. They mainly originate from the evaluation of the sublead- 
ing contribution to the vertex function and uncertainties in excited meson wave functions. 
In Tables IVllVIII we compare our results for form factors of the weak Bs decays to scalar 
Kq, axial vector Ki and tensor K2 mesons at maximum recoil point g^ = with previous 



theoretical calculations |l7l. |22|-|29||. Predictions are mostly available for the weak Bs decays 



to scalar K^ mesons. Such transitions were studied in the light-cone sum rules [22|, |2J] , QCD 



sum rules [23|, |26| and perturbative QCD approach [25|, |27|. We find that our result for the 



form factor r+(0) is consistent with Ref. [23|, but lower than in other calculations | 22l . 2j 
Contrary we predict larger absolute value for the form factor r„ (0) than in Refs. J22|, |24 
Predictions for form factors of the weak Bg decays to axial vector Ki mesons are compared 
in Table IVTl The light-cone sum rule approach is used in Ref. |28|, while Ref. [l7| employs 



27|. 
26|. 



perturbative QCD. From this table we find that our and these theoretical approaches predict 



11 



TABLE V: Comparison of theoretical predictions for the form factors of the weak Bg — )• Kq 
transitions at maximum recoil point q'^ = 0. 
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TABLE VI: Same as in Table |V] but for the Bs -^ Ki transitions. 
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5a(0) -0.33 ±0.03 -0.15l[!:[J? 0.03 ± 0.01 

5^1 (0) -0.08 ±0.01 -0.1ll[j:[|l 0.11 ±0.07 

/ia(0) 0.29 ±0.03 0.49 ±0.08 0.20 ± 0.05 

/iVi(O) 0.08 ±0.01 0.38 ±0.06 0.87 ±0.25 



significantly different values of form factors at maximum recoil. One of the origins can be 

the adopted values of the mixing angle </?, defined in Eq. fl28|) . of axial vector Ki mesons. 

Note that in our model this angle yj is explicitly calculated by diagonalizing the mass matrix 

5|, while in Refs. IT], |28| different phenomenologically motivated values are used. The 



form factors of the weak Bg decays to tensor K2 mesons were calculated in Ref. 29| in the 
perturbative QCD approach. The obtained values of the form factors ty(0) and tyii(O) are 
slightly lower than in our model. 

V. CHARMLESS SEMILEPTONIC B^ DECAYS 

We can now apply the calculated form factors for the evaluation of the semileptonic Bg 
decays to ground state and orbitally excited K mesons. The differential decay rate of the 
Bs meson to a, K [K^*' or Kj) meson can be expressed in the following form [3C 



K 



ub\ 



dq" (27r)3' ""' 24M|g 



2 



(29) 
where Gp is the Fermi constant, V^b is the CKM matrix element, A = A(M^,M|.,g^) = 
M% + M^ + q^ - 2(M^M| ± M^q^ + M^q"^), mi is the lepton mass. The transverse H±, 
longitudinal Hq and time Ht helicity components of the hadronic tensor are defined through 
the transition form factors calculated in the previous sections. The corresponding relations 
for decays to ground and orbitally excited mesons are given in Appendices C and E of 
Ref. |3|, respectively. 





TABLE VII: Same as in Table |V] but for the Bg -^ K^ transition. 




F(0) 


this paper 


[29] 


V(0) 

tAM 


-0.34 ±0.03 
-0.17 ±0.02 


-o.i8l°:°l 
-0.111°:°^ 
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FIG. 5: Predictions for the differential decay rates of the semileptonic Bs — t- K^*'lvi decays. 



TABLE VIII: Comparison of theoretical predictions for the branching fractions of semileptonic 
decays Bs -^ K^*Hvi (in 10"^). 



Decay 



this paper 



[21] 



1.64 ±0.17 
0.96 ±0.10 
3.47 ±0.35 
1.67 ±0.17 



[19] 



Bs 
Bs 
Bs 
Bs 



> Keve 



1.27 
0.778 



0.49 
-0.30 

0.268 
-0.201 



1.47 ±0.15 

1.02 ±0.11 
2.91 ±0.26 
1.58 ±0.13 



A. Semileptonic Bs decays to ground state K^ "i mesons 

First we calculate the rates of the semileptonic Bg decays to ground state K^*^ mesons. 
Substituting the form factors obtained in Sec. [TTll in Eq. (!29|) we get the corresponding dif- 
ferential decay rates. They are plotted in Fig. [5] both for decays involving electron e and 
r lepton. Integrating these differential decay rates over q^ we find the total decay rates. 
In calculations we use the value of the CKM matrix element \V,i.h\ =_(4.05 ± 0.20) x 10~^ 



found previously from the comparison of the predictions of our model [ij, [SlJ with measured 
semileptonic B — )■ T[{p)lvi decay rates. The kinematical range accessible in the heavy-to- 
light Bs — )■ K^*'^ transitions is very broad, that makes the knowledge of the q^ dependence 
of the form factors to be an important issue. Therefore, the explicit determination of the 
momentum dependence of the weak decay form factors in the whole q^ range without any 
additional assumptions is an important advantage of our model. The calculated branching 
fractions of the semileptonic Bs -^ K^*Hui decays are presented in Table rvTTTl in comparison 
with other theoretical predictions 19|, |2l|. The perturbative QCD factorization approach is 



used in Ref. [2l|, while in Ref. [19| light cone sum rules are employed. From the compari- 
son in Table IVIIII we see that all theoretical predictions for the Bs semileptonic branching 
fractions agree within uncertainties. This is not surprising since these significantly different 
approaches predict close values of the corresponding weak form factors (see Table [In|) . 

We can use the calculated values of the semileptonic Bs decay branching fractions to 
obtain predictions for the ratios of such decay involving r lepton and electron or muon: 
R{K) = Br{Bs -> KTUr)/Br{Bs -> KeUe) = 0.59 ± 0.05 and R{K*) = Br{Bs -^ 
K*tVt-)/ Br{Bs — )■ K*eVf,) = 0.48 ±0.04. The interest to such ratios is stimulated by recently 
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TABLE IX: Comparison of theoretical predictions for the branching fractions of semileptonic decays 



Bs 


— > 


Kyiui (in 10^ 


'% 


















Decay 


this paper 


[23] 


[24] 


[25] 


[28] 


[17] 


[29] 




Bs 


— > KQCUf. 


0.71 ±0.14 


o.36l°:i 


i-3tJi 


2.45tl:^^ 










Bs 


-^ K^TVr 


0.21 ±0.04 




0.52l°:?3^ 


1.09l°i? 








Bs 


-^ 


Ki{1270)eiye 


1.41 ±0.28 








4.531^:^^ 


"-••'^-2.89 




Bs 


—?- 


Ki{1270)TUr 


0.30 ±0.06 










2mt\f, 




Bs 


— ^ 


Ki{UOO)ei^e 


0.97 ±0.20 








^mt\^l 


0.031°:°^ 




Bs 


Bs 
Bs 


Ki(1400)ri/^ 
— ;> K2ei'e 

-^ K2TUr 


0.25 ±0.05 
1.33 ±0.27 
0.36 ± 0.07 










o.oitn? 


0.73l°:i 
0.25t°:}^ 



found deviations of experimental data from theoretical predictions for the similar ratios for 
the semileptonic B decays to D mesons R{D^*^) = Br{B — )■ D^*Wh'r)/Br{B — )■ D''*^eve) 
(see, e.g., discussion in [31] and references therein). 

Summing up different contributions listed in Table IVIIII we obtain the prediction for 
the total semileptonic Bs decay branching fraction to the ground state K mesons to be 
Br{Bs -^ K^*hve) = (5.11 ± 0.51) x 10"^ and Br{Bs -^ K^*hur) = (2.63 ± 0.26) x 10"^ 

B. Semileptonic Bg decays to orbitally excited Kj mesons 

Now we calculate the branching fractions of the semileptonic Bs decays to orbitally excited 
Kj mesons. We substitute the form factors obtained in Sec. IIVI in the expression for the 
differential decay rate fl29|) . The resulting decay rates are plotted in Fig. [61 Integration 
over g^ gives the total semileptonic Bs — ?■ Kjlvi branching fractions which are given in 
Table IIXI We see that our model predicts close values (about 1 x 10~^) for all semileptonic 
Bs branching fractions to the first orbitally excited Kj mesons. Indeed, the difference 
between branching fractions is less than a factor of 2. This result is in contradiction to the 
dominance of specific modes (by more than a factor of 4) in the heavy-to-heavy semileptonic 
B — )• Djlui and Bs — )■ D^Jjlui decays |3|, |32|, but it is consistent with predictions for the 



corresponding heavy-to-light semileptonic B decays to orbitally excited light mesons [33 
The above mentioned suppression of some heavy-to-heavy decay channels to orbitally excited 
heavy mesons was mostly pronounced in the heavy quark limit and then slightly reduced 
by the heavy quark mass corrections which are found to be large [32J] . Thus our result once 
again indicates that the s quark cannot be treated as a heavy one and should be considered 
to be light instead, as we always did in our calculations. 

In Table HXl we compare our predictions for the semileptonic Bs branching fractions to or- 
bitally excited Kf mesons with previous calculations [l7|, |23|-|25|. l28l . |29| . The consideration 



in Ref. [23[ is based on QCD sum rules. The light cone sum rules are used in Refs. |24l . I28lj . 
while Refs. IT], |25|, |29| employ the perturbative QCD approach. Reasonable agreement be- 
tween our results and other predictions |23|, |2J, |29| is observed for the semileptonic Bs decays 



to the scalar and tensor K mesons. The values of Ref. J25| are almost a factor 3 higher. 
For the semileptonic Bs decays to axial vector K mesons predictions are significantly dif- 
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FIG. 6: Predictions for the differential decay rates of the semileptonic B — t- Kj 'Ivi decays. 



ferent even within rather large errors. Therefore experimental measurement of these decay 
branching fractions can help to discriminate between theoretical approaches. 

For the total semileptonic Bg decays to first orbital excitations of K mesons we get 
Br{Bs -^ K^j^ePe) = (4.4 ± 0.9) x 10"^ and Br{B, -^ K^j^rUr) = (1-1 ± 0.2) x 10"^. These 
values are close to the ones found for the semileptonic Bg decays to ground state K mesons 
in Sec. IV A[ The similar pattern of the branching fraction dependence on the excitation 
of a final meson was previously found for the heavy-to-light semileptonic B decays |33| . 
On the other hand, for the heavy-to-heavy semileptonic B ^ D and Bg — )■ Dg decays the 
pronounced hierarchy, where the decay branching fractions rapidly decrease with the growing 
excitation of the final heavy {D or Dg) meson, is observed [3|, |32 . 



VI. CHARMLESS NONLEPTONIC DECAYS OF B, MESONS 



We further apply calculated Bg decays form factors to the evaluation of the charmless 
nonleptonic decays of the Bg meson. In the discussion we closely follow our previous con- 
sideration of nonleptonic decays of B mesons to ground and orbitally excited states of light 
mesons in the factorization approximation 33|, |3J]. To simplify the problem we limit our 
analysis to the calculation of the decay processes dominated by the tree diagrams. For such 
decays the matrix elements of the effective weak Hamiltonian Hes, governing nonleptonic 
decays Bg — > i^~'°M"'"'°, where M is a light (vr or p) meson, can be approximated by the 
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product of one-particle transition amplitudes 

Gr p 



{K-M'^lH.sm 
{K^M'\H,s\Bs) 



^V:,Vudaf{K-\{bu)v-A\Bs){M+\{du)v-A\0) 



(30) 



with 



,eff 



Oi 



[a4 + aio + rg{aQ + as)], 



af 



a2 — 



VtlVt, 



3 3 1 

-04 ip -a-j + -ag + -aio 



r,(2a6 - ag) 



v:,v^d 

where terms in square brackets result from the penguin contributions, which are small nu- 
merically; minus and plus correspond to M° = vr" and M° = p°, respectively. The quantities 
a2n-i = C2n-i + C2n/^c and 02™ = C2„ + C2„_i/A^c iji = 1,2... and A^c is the number of colors) 
are combinations of the Wilson coefficients Cj, which we take from Ref. [35|, and r^ can be 
found, e.g., in Ref. j36l] . 

The matrix element of the weak current J^ between vacuum and a final pseudoscalar 
(P) or vector (V^) meson can be parametrized by the decay constants fpy 

{P\qil^im2\Q) = ^fpp'i>, (^|gi7M?2|0) = e^Mvfv (31) 

As a result the corresponding nonleptonic matrix element factorizes in the product of the 
weak Bs ^^ K decay form factors, which were calculated in the previous sections, and 
decay constants. The pseudoscalar fp and vector fy decay constants of light and heavy 
mesons were calculated within our model in Ref. |37|. Their values are in agreement with 
the available experimental data [2| . For the calculations we use the following values of the 
decay constants: /^r = 0.131 GeV, fp = 0.208 GeV. The relevant CKM matrix elements are 
iKdl = 0.975, \Vtd\ = 0.0087, \Vti,\ = 0.999 0. 

In Table |X] we present the predictions for the branching fractions of the charmless non- 
leptonic Bg decays to ground state K mesons obtained in the factorization approximation 
with our model form factors. There we also give results of other theoretical approaches 
|l6l . I20I |38|-|40| and available experimental data J2|. Calculations in Ref. |l6| are done in 
perturbative QCD approach. QCD factorization is employed in Refs. 38|, |40|. The authors 
of Ref. |39| use soft-coUinear effective theory, while considerations in Ref. 20| are based 
on an approximate six-quark effective Hamiltonian. Experimentally only the Bg — )■ K^n^ 
branching fraction was measured |2| and for the B^ — )■ K*^p^ decay upper limit is available. 
All theoretical predictions agree with data. In fact all theoretical results for decays into 
charged K^*^~ and 7r+(p"'") mesons are consistent within rather large error bars, while for 
decays involving neutral mesons deviations are larger. 

We apply the same factorization approach for the calculation of the nonleptonic decays 
to orbitally excited Kf mesons. Using form factors obtained in Sec. IIVI we get predictions 
for the nonleptonic branching fractions and present them in Table IXII in comparison with 
previous estimates 27|, |4l|. A few predictions are available for selected modes. In Ref. |27| 
decay Bg — )■ Kq~ p'^ was considered within the perturbative QCD factorization approach. 
The obtained central value of this decay branching fraction is almost a factor of 4 larger 
than our result. This is the consequence of a 2 times larger form factor r+(0) in Ref. J27| 
than in our model (see Table |V|) . Decays involving the tensor K2~ meson were considered in 
Ref. [4l| within the Isgur-Scora-Grinstein-Wise II model. The predictions are approximately 
a factor of 2 lower than our central values of branching fractions. 



TABLE X: Comparison of 
tonic Bg decays to ground 
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various predictions for the branching fractions of the charmless nonlep- 
state K mesons with experiment (in 10^^). 





Decay 


this paper 


[16] 


[38] 


[39] 


[40] 


[20] 


Experiment [2] 


Bs 


— )• K"7r+ 


8.7 ±2.7 


7.6+0 


10.2+fO 


4.9 ±1.8 


c; q+0.5 


7.lt?i 


5.3 ±1.0 


Bs 


^ K~p+ 


24.0 ± 7.2 


I7.8tli 


24.5l:lli 


10.2 ±1.0 


i4.7t^:3^ 


17.6l|2 




Bs 


-^ K*-TT+ 


8.6 ±2.6 


7-6i1i 


8.7tli 


6.6 ±0.7 


7.8ig| 


7 9+5.6 
'•^-2.3 




Bs 


-^ K*-p+ 


25.4 ±7.6 


20.9t^:^ 


25.2tti 




21.61^2 


2i.o+^;'5' 




Bs 


^K^TT^ 


0.25 ±0.08 


16+0-11 
U-lO-0.06 


o.49+°:i 


0.76 ±0.41 


1-71L^ 


1 1+°-'^ 




Bs 


: -^ K^pO 


0.67 ±0.20 


n no+0.07 


n CI +1.26 

"•o-i-o.eo 


0.81 ±0.09 


i.9l:i:i 


0.610:3 




Bs 


^ ET^OvrO 


0.24 ±0.07 


o.o7l°:°i 


n 9^+0-46 


1.07 ±0.19 


0.89tlil 


n ■^+0-2 




Bs 


^ K*Op^ 


0.71 ±0.21 


0.33l°:ll 


1 r:+3.3 




1.3lf, 


1 n+0-4 

l-"-0.3 


<767 



TABLE XI: Branching fractions of the nonleptonic Bs decays to orbitally excited Kj mesons (in 

10-6). 



Decay 



this paper 



[27] 



[41] 



Bs 



K*-v^^ 



Bs ^ K*-p+ 
Bs ^-firi(1270)7r+ 
Bs^Ki{l270)p+ 
Bs-^ K{{UOO)tt^ 
Bs ^Kf (1400)p^ 

Bs -^ k;^tt+ 



Bs 



K;-p^ 



9.6 ±3.8 
27 ±10 
29 ±12 
76 ±30 

21 ±8 

54 ±21 

17 ±6 

47 ±18 



108 



-34 
-31 



7.8 
23 



VII. CONCLUSIONS 



The form factors of the Bs weak decays to strange mesons were calculated in the frame- 
work of the relativistic quark model based on the quasipotential approach. Decays both to 
the ground K^*'> and to orbitally excited Kj mesons were considered. The form factors were 
determined as the overlap integrals of the related meson wave functions, found in the previ- 
ous meson mass spectrum calculations [4, [Sj, in the whole broad kinematical range without 
additional assumptions about their g^ dependence. The important relativistic effects, such 
as transformations of the meson wave functions from the rest to a moving reference frame 
and contributions of the intermediate negative-energy states, were consistently taken into 
account. 

We used these form factors for the evaluation of the charmless semileptonic Bs decay rates 
to ground state and orbitally excited K mesons. It was found that total branching fractions 
of semileptonic Bg decays to ground and first orbitally excited K mesons have close values 
about 5x 10~^. Summing up these contributions we get (9.5±1.0) x 10~^. This value is almost 
2 orders of magnitude lower than our prediction for the corresponding sum of branching 
fractions of the semileptonic Bg to Dg mesons [Sj as it was expected from the ratio of CKM 
matrix elements \Vub\ and | V^b|. Therefore the total semileptonic Bg decays branching fraction 
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is dominated by the decays to Dg mesons and in our model is equal to (10.3 ± 1.0)% in good 
agreement with the experimental value Br{Bs — )■ Xez/e)Exp. = (9.5±2.7)% [2|]. The obtained 
predictions for charmless semileptonic Bg decays were compared with previous calculations. 
It was found that different theoretical approaches yield close values of branching fractions 
for decays to ground state K^ "^ mesons agreeing within uncertainties, while the ones for 
the decays to orbitally excited Kj mesons differ significantly from each other. The latter 
observation can help to discriminate between theoretical models. 

Charmless two-body nonleptonic Bg decays dominated by tree diagrams were considered. 
The factorization approximation was used which allowed to express the decay matrix ele- 
ments as the product of weak transition matrix elements and decay constants. The branching 
fractions of the nonleptonic decays to ground state K^*'' or orbitally excited Kj meson and 
pion or p meson, both charged and neutral, were calculated. The obtained results were 
confronted with previous theoretical predictions and experimental data, which are available 
for only few of the considered decays. Good agreement with data and other evaluations is 
found for the nonleptonic decays to the ground state K^*' meson and a light meson, while 
again for decays involving the orbitally excited Kj mesons significant disagreement between 
predictions of different approaches is observed. 
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